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Abstract. We introduce a new notion of the determinant, called symmetrized de- 
terminant, for a square matrix with the entries in an associative algebra A. The 
monomial expansion of the symmetrized determinant is obtained from the standard 
expansion of the commutative determinant by averaging the products of entries of 
the matrix in all possible orders. We show that for any fixed finite-dimensional asso- 
ciative algebra A, the symmetrized determinant of an n x n matrix with the entries 
in A can be computed in polynomial in n time (the degree of the polynomial is linear 
in the dimension of .4). Then, for every associative algebra A endowed with a scalar 
product and unbiased probability measure, we construct a randomized polynomial 
time algorithm to estimate the permanent of non-negative matrices. We conjecture 
that if A = Mat(d, M) is the algebra of d x d real matrices endowed with the standard 
scalar product and Gaussian measure, the algorithm approximates the permanent of 

a non-negative n x n matrix within 0(7?) factor, where lim^ >+OQ 7^ = 1. Finally, 

we provide some informal arguments why the conjecture might be true. 



1. Introduction 

(1.1) Permanent and determinant. Let A = (a^) be an n x n real matrix and 
let S n denote the symmetric group of all permutations a of the set {1, . . . , n}. The 
number 

n 

per A = ^ Yl a ia(i) 

cr£S n i=l 

is called the permanent of A. We will be interested in the case when A is non- 
negative: ciij > for i, j = 1, . . . , n. 

To compute or to approximate the permanent efficiently is one of the most 
intriguing problems, see, for example, [Jerrum and Sinclair 89], [Karmarkar et al. 
93], Chapter 18 of [Papadimitriou 94] , [Linial et al. 98], [Barvinok 97] and [Barvinok 
99]. 
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On the other hand, the determinant 



n 

detA= ( s S n a ) II a M*) 

a£S n i=l 

can be computed using 0(n 3 ) arithmetic operations (even when A is a matrix over 
a commutative ring), see, for example, Sections 11.1 and 15.1 of [Papadimitriou 
94]. 

Hence it seems natural to try to use determinants to approximate the permanent. 

(1.2) Permanent estimators. The following approach is due to Godsil and Gut- 
man and independently to Girko (see, for example, Chapter 8 of [Lovasz and Plum- 
mer 86] and Section 3 of Chapter 2 of [Girko 90]. Let us choose a probability 
distribution |ionl with the properties that 

E x = / x d/j,(x) = and Ei 2 = / x 2 dfi(x) = 1. 
Jr Jr 

Let us sample n 2 numbers My : i,j = 1, . . . , n independently and at random from 
fx. Let B = {pij) be the matrix defined by 

b{j — U{j (Xij 

and let 

a = (detB) 2 . 

Then a is a random variable and it turns out that E a = per A. Hence, in principle, 
per A can be approximated by averaging sufficiently many determinants. It turns 
out that by extending the ground field, one can ensure a better concentration 
of a around its expectation and hence a better computational complexity of the 
approximation. 

In [Karmarkar et al. 93] it was shown that if one allows My to be complex 
numbers (we must define then a = |detS| 2 ), one can make the variance of a 
smaller. More precisely, the variance of a for the distribution // that chooses the 
cubic roots of unity with the probability 1/3 each is exponentially smaller in the 
worst case than that for \i that chooses —1 and 1 with the probability 1/2 each. 
In [Barvinok 99], it was shown that if fi is a Gaussian distribution in R then 
with high probability a random value of a approximates its expectation within a 
c n factor for c ~ 0.28; similar behavior is observed when fi is a complex Gaussian 
distribution, but in the complex case the constant c gets better: c w 0.56. Moreover, 
if n is a quaternionic Gaussian distribution (in which case a should be the Study 
determinant of £?, that is the determinant of the complexification of £?, see, for 
example, [Aslaksen 96]), the constant c gets even better: c ~ 0.76. 

This suggests that it may be of interest to construct more determinant-type 
estimators, as they could provide still better concentration properties of a. 
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Let us summarize some useful properties of the estimator. The estimator a = 
a(uij) : i,j = 1, . . . , n is a function of n 2 independent and identically distributed 
random variables Uij, such that: 

(1.2.1) For each choice of Uij, the value of a(uij) can be computed using a polyno- 
mial in n number of arithmetic operations; 

(1.2.2) The expected value of a is per A; 

(1.2.3) For any choice of (uy), the value of a(uij) is a non-negative real number; 

(1.2.4) Let us fix all Uij except those in one row i$ (resp. in one column jo). Then 
a is a quadratic form in {ui j : j = 1, . . . , n} (resp. in {uij Q : i = 1, . . . , n}). 

In this paper, for every associative finite-dimensional algebra A we construct an 
estimator with the properties (1.2.1)-(1.2.4). 

(1.3) Non-commutative determinants. Exploring analogues of determinants 
over non-commutative rings and algebras is a very old topic, see, for example, 
[Aslaksen 96] for a survey and [Gelfand and Retakh 97] for new developments. In- 
terestingly, most important non-commutative determinants, such as the Dieudonne 
determinant, quasideterminants of Gelfand and Retakh and the Moore determi- 
nant of a Hermitian quaternionic matrix turn out to be computationally efficient. 
However, with the notable exception of the Moore determinant, they are not close 
enough to the monomial expansion of the commutative determinant to produce a 
permanent estimator. On the other hand, the most straightforward version due to 
Cay ley 

CdetA= ^ ( s g ncr ) a io-(i) • • -a-nain) 

for an n x n matrix A = (a^) with the entries a^- in an associative algebra A 
would have suited our purpose perfectly well had there been any reason to believe 
that this expression can be efficiently (in polynomial in n time) computed, see also 
Section 5.2. 

The ambiguity that prevents the straightforward extension of the determinant 
to non-commutative algebras is that there is no natural choice for the order of the 
factors in the product Y\7=i a io{i)- We resolve this ambiguity by taking the average 
of the products in all possible n\ orders. Hence for annxn matrix A = (a^) with 
the entries in an associative algebra A over a field F of characteristic 0, we write 

Sdet A = — } ^ ( S S n a ) ( S S n T ) a v(l)T(l) ■ ■ ■ a a (n)r(n) 
(cr,r)eS„ XS„ 

("sdet" stands for the "symmetrized determinant"). It is easy to see that if A 
is commutative, we get the standard determinant. We prove that for any fixed 
finite-dimensional algebra A, the value of sdet A for an n x n matrix over A can be 
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computed in a polynomial in n time. More precisely, if the dimension of A as an 
F-vector space is r, sdet A can be computed in 0(n r+3 ) time. 

The paper is organized as follows. In Section 2, we discuss mixed discriminants, 
which are crucial for our proof in Section 3 of the polynomial time computability 
of the symmetrized determinant. In Section 4, for every finite-dimensional associa- 
tive algebra A over R endowed with a scalar product and an unbiased probability 
distribution [i, we construct a permanent estimator which satisfies (1.2.1)-(1.2.4). 
In Section 4, we conjecture that if A = Mat(<i, M.) is the algebra of d x d matrices 
endowed with the standard scalar product and Gaussian probability measure, then 
the algorithm approximates the permanent of an n x n matrix within a 0(^2) f ac ~ 

tor, where lim^ >+00 ^y d = 1. We also provide some intuitive argument supporting 

the conjecture. 

2. Preliminaries: Mixed Discriminants 

(2.1) Definition. Let Ai, . . . , A n benxn matrices over a field F of characteristic 
0. We write Ak = (a^j), where a^- G F for i,j = 1, . . . ,n and k = 1, . . . ,n (k 
is the index, not the power). Let t±, . . . ,t n G F be variables. The expression 
det(ti^4i + . . . + t n Anj is a homogeneous polynomial in t±, . . . , t n of degree n and 
its normalized coefficient 

1 d n 

D(A 1 , ... ,A n ) = — — _ det(Mi + . . . + t n A n ) 

is called the mixed discriminant of Ai, . . . , A n . In terms of the entries (a^) of the 
matrices Ai, . . . , A n , the mixed discriminant can be written as 



1 n 
(2.1.1) D(A l ,...,A n ) = - (sgna)(sgnr)n«" 

' (a,T)eS n xS n k=l 



Mixed discriminants are symmetric, that is, 

(2.2) D(A U ...,A n ) = D(A H1) , . . . , A Hn) ) 

for any permutation : {1, . . . , n} — > {1, . . . , n}. For various properties of mixed 
discriminants, see, for example, Section 5.2 of [Bapat and Raghavan 97]. 

We are particularly interested in the situation when the number of different ma- 
trices among A\, . . . , A n is small. The following result is Lemma 9.3 from [Barvinok 
97]. For the sake of completeness, we present its proof here. 

(2.3) Lemma. Let ki, . . . ,k r be non-negative integers such that k\ + . . . + k r = n 
and let Ai, . . . , A r be n x n matrices. Then 

D [Ai, . . . , A\, . . . , A r , ... , A r 




0<m r <k r 




£ +m '(™ 1 1 ) "■(£) ^A 1 + ... + mr A T ). 



Proof. For a subset u C {1, . . . , n}, let 

1 if % G a; 



if % oo 



and let t w = (ti(cu), . . . , t n (a>)) be the indicator of a;. One can observe that if p is 
a homogeneous polynomial of degree n > in n variables £i, . . . , £ n , then 

^L_ P =(_!)• £ (-l)'"'p(t„). 

wC{l,... ,n} 

Indeed, it suffices to check the identity for monomials t" 1 ■ ■ If some oti = 

then the right hand side is since the terms corresponding to w U {i} and u \ {i} 
annihilate each other. For the monomial t\ • ■ -t n the right hand side is 1 since the 
only non-zero term is for u = {1, . . . , n}. 

Given n x n matrices A±, . . . , A n , let us apply the above identity to 

p(t 1 , ... ,t n ) = det(tiAi + . . . + t n A n ). 

Suppose that the set {Ai, . . . , A n } consists of k\ copies of A\, . . . , k r copies of A r , 
where A\, . . . ,A r are distinct, and let Si = {j : Aj is a copy of A^ for i = 1, . . . , r. 
Then \Si\ = ki, Si U . . . U S r is a partition of the set {1, . . . , n} and 

p(t w ) = det(mi^4i + . . . + m r A r ), where rrii = \ui D Si\. 

Moreover, for given mi, . . . , m r , there are exactly ( 1 J • • ■ ( r ) subsets u with 

\m 1 J \m r J 

rrii = \uC\Si\, since for each i = 1, . . . , r, we have to choose rrii elements of u> from 
Si for alH = 1, . . . , r independently. The proof now follows. □ 

(2.4) Corollary. Suppose that r is fixed. Given n x n matrices Ai, . . . , A r , com- 
puting 

D ( ^4-i, ■ • • , A\j . . . , A r , . . . , A r 

ki copies k r copies 

using the formula of Lemma 2.3, takes 0(n r+3 ) arithmetic operations. 

Proof. The number of summands does not exceed (n + l) r (a better estimate is 
(n/r + l) r ) and the determinant of an n x n matrix can be computed using (9(n 3 ) 
arithmetic operations. □ 

3. Symmetrized Determinant of a Matrix over an Algebra 

Let A be an associative algebra over F, where F is a field of characteristic 0. 
Hence A is a vector space over F with an addition "+" and associative (but not 
necessarily commutative) multiplication •. For example, one can choose A to be 
the algebra Mat(<i, F) of all d x d matrices over F. We will assume that as a vector 
space, A is finite-dimensional, dimp^A < oo. 
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(3.1) Definition. Let A = (a^), a^- G A be an n x n matrix over A We call 

Sdet A = — } ^ ( S S n CT ) ( S S n r ) a <T(l)r(l)a (J (2)r(2) • • • Oa(n)T(n) 

(ff,r)6S„xS„ 

the symmetrized determinant of A. In other words, sdet A is obtained by taking a 
diagonal of A (that is, picking one entry from each row and column of A), multi- 
plying the entries on the diagonal in all possible orders, taking the average of the 
resulting n\ products and adding that average with the appropriate sign found by 
the same rule as for the usual commutative determinant. Hence sdet A G A. We 
denote by Mat(n, .4.) the set ofnxn matrices A = (a^ ) with a^- G A. 

Remark. If A is commutative, the expressions 

~. Yl ( S g n(T )( S g nr ) a f T(l)r(l)---a f T(n)r(n) 

(<r,T)es n xS n 

and 

coincide but if A is not commutative, they may differ. 

(3.2) Theorem. Let e\, . . . ,e r span A as an ¥ -vector space. Let A = (a^) be an 
n x n matrix over A, so 

A = A\e\ + . . . + A r e r , 

where Ai, . . . ,A r are n x n matrices with the entries in ¥. 

For an r -tuple k\, . . . , k r of non-negative integers such that k± + . . . + k r = n, 
let <&(fci, . . . , k r ) be the set of all maps 4> : {1, . . . , n} — ► {1, . . . , r} such that 
|(/> -1 (l)| = ki, . . . , |(/> _1 (r)| = k r . Let us define elements u(ki, . . . , k r ) G A by 

u{k u . . . , k r ) = Y e <Ki) e <K2) - ' - e 4>{n)- 
<£e«i>(/ci,... ,k r ) 

Then 

sdet A = ^ D(Ai, ... ,Ai,... , A r , . . . , A r )u(ki, ... , k r ). 



k~i_ j - - - •jfc'T' ^ 
ki+...+k r =n 



ki copies k r copies 



Proof. Let A^ = (a^), where af- G F for i, j = 1, . . . , n and k = 1, . . . , r. Hence 

Er 
k=l a ij e k- 
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By Definition 3.1, 

sdetA=^ (sg n ^)(sg n ^)a<T(i)r(i)a fT (2)r(2)---a fT ( n)T ( n ) 

(a,T)es n xs„ 

^ r r 

= ~\ Z (sgna)( S gnr)(j]a^ 1)r{1) e,) ... (XX(»Mn) e *) 

' (<j,r)65„x5„ fc=l fc=l 

Let $ be the set of all maps : {1, . . . , n} — > {1, . . . , r}. Then 

r r 

k \ f\~^ k \ \^ </>(l) <^>(n) 

(,2^ a <r(l)r(l) e fc J'-' {2_^ a <7(n)r(n) e k ) ~ 2^ a <r(l)T(l) ' ' ' a *(n)r(n) e Hl) ' ' ' e <M™) ■ 

k=i k=i 4>e^> 

Hence 

SdetA = ^! ^ (sgna)(sgnr) ^ a^ r(1 y ■ ■ a^^e^ ■ ■ ■ e Hn) 

' (a,r)eS n xS n 06* 
= 5Z^ ( S g na )( S g nr ) a ^!)r(l)"-- a ^n)r(n)^(l)---^ 

^e* ' (a,T)eS n xS n 
= Z -°(^(i)' • • • ' AKn))e<Ki) ' ' ' e <£(n)- 

Now, set $ is a disjoint union of the sets $>(k\, . . . , k r ) and by the symmetry of the 
mixed discriminant (see (2.2)), for any <fi e 3>(fci, . . . , fc r ) we have 

-D(^(i), • • • A/>(n)) = D( Ai, . . . , Ai , . . . , A r , . . . , A r y 

fci copies k r copies 

Summarizing, we get 

sdetA= Y D ( Au ■ • • , •■■, Av • • ; A r ) ^ e^y ■ e^ n) 

£+"+k^!n *i copies k r copies 0e*(fci,... ,fc r ) 

and the proof follows. □ 

Next, we present an algorithm for computing the symmetrized determinant of a 
matrix over A. We assume that there is a basis ei, . . . , e r of A as a vector space 
and that the entries of A G Mat(n, A) are given by their coefficients in the basis: 
ciij = ajjCi + . . . + aL-e r . The multiplication in A is described by the structural 
constants that are the scalars from F such that 



r 



(3.3) e i e j = Y,Pfak- 



fe=i 



We need a simple Lemma. 



(3.4) Lemma. Let u(k\, . . . , k r ) G A be the elements defined in Theorem 3.2. Let 
I = {i:ki> 0}. Then 

u(ki, . . . , k r ) = ^ ^ u{ki, • • • , ki—i, ki — 1, . . . , k r )ei. 
iei 

Proof. For G <&(fci, . . . , fc r ) we have 

e^(i) • • -e^( n ) = ^ ( i) • • -e^( n _i)^e^( n ), where 0(n) G I. 

If 0(n) = z then the restriction <f>' : {1, . . . , n — 1} — > {1, . . . , r} belongs to the 
set . . . , fcj — 1, ki+i, . . . , fc r ). Vice versa, every <fi' G . . . , fcj — 

1, /ci+i, . . . , fc r ) extends to e <&(fci, . . . , fc r ) by letting 0(n) = z. The proof now 
follows. □ 

(3.5) Algorithm for computing sdet A for A G Mat(n, .4) 

Input: An algebra .4 given by the structural constants {/&} in some F-basis 
ei, . . . ,e r of A and an n x n matrix A G Mat(n,*4) given by n x n matrices 
4.1, . . . , 4. r G Mat(n, F) such that A = A x e x + . . . + A r e r . 

Output: The element sdet A G A given by its coefficients in the basis ei, . . . , e r . 

Algorithm: First, for all r-tuples of non-negative integers k±, . . . , /c r such that + 
. . . + k r = n, using the identity of Lemma 2.3, we compute the mixed discriminant 

D(Ai, . . . , Ai, . . . ,A r , . . . , A 7 

fci copies fc r copies 

Now, using Lemma 3.4 recursively, we compute u(ki, . . . , k r ). We start with 

u(0,...,0, 1, 0,...,0) = e, 

i-th position 

and applying Lemma 3.4 and (3.3), successively compute the expansions of 
u(ki, . . . , k r ) with ki + . . . + k r = 1, . . . , n in the basis e±, . . . ,e r . Finally, we 
use Theorem 3.2 to compute sdet A. 

An important observation is that if dimjr.4 = r is fixed, the algorithm has 
polynomial time complexity. 

(3.6) Theorem. For a fixed r, given an n x n matrix A G Mat(^4, n), Algorithm 
3.5 computes sdet A G A using 0(n r+3 ) arithmetic operations. 

Proof. Theorem 3.2 and Lemma 3.4 imply that the algorithm indeed returns the 
correct value. The number of all r-tuples of non- negative integers k\ , . . . ,k r such 
that ki + . . . + k r < n is ( n ^ r ) , which is a polynomial in n of degree r. The 
complexity bound follows from this and Corollary 2.4. □ 



4. Permanent Estimators 



Let us fix an IR-algebra A. We assume that as a real vector space, A is finite- 
dimensional, so as a vector space A is isomorphic to M. m for some m. Suppose that 
there is a scalar product (•, •) : A x A — > ffi. In other words, for every two elements 
a, b G A, a real number (a, b) is defined, such that 

(a, b) = (b, a); 

(aa + /?&, c) = ai(a, c) + c) for a,b,c E A and ct, /3 G K. and 
(a, a) > for all a E A. 

For example, if A = Mat(<i, R) is a matrix algebra, one may choose (a, 6) = 
T^afr*). Generally, we don't assume any relation between the multiplication in A 
and the scalar product. We also allow the form (•, •) to be degenerate, that is, we 
allow (a, a) = for some non-zero a. As usual, we define ||a|| 2 = (a, a). 

Suppose further, that there is a Borel probability measure fx on A with the 
properties: 




and for any (non-commutative) polynomial 

p(x!,... ,x n ) = lii,...,i n x h'-- x i n i where 7ii,...,i„eR 

l<il,... ,i n <n 

we have 

(4.2) E (p(xi,... ,x n ),p(xi,... ,x n )) < +oo, 

provided xi, . . . , x n E A are chosen independently and at random. 

(4.3) Theorem. Let A be an M,-algebra with a scalar product (■, ■) and a probability 
measure p satisfying (4-l)-(4-%)- Let us define a constant c{n,p) > as follows: 

c{n,p) = E ||Z|| 2 , where Z = — ^ Y a{1) ■ ■ ■ Y a{n) 

and Yi, . . . ,Y n are sampled independently and at random from A. 

For a given real non-negative matrix A = (dij), let us define a random n x n 
matrix B e Mat(n, A), B = (bij) as follows: 



where X^ are sampled independently and at random from A. Then 

E || sdetS|| 2 = c(n, fj,) per A. 

Proof. By Definition 3.1, 



1 xn .1/2 

sdetS = — ^ (sgna)(sgnr) f JJa (T ( i ) T ( i ) J X a{1)r{1) ■ ■ ■ X (j(n)r(n) . 



1/2 

>gna)(sgnr)( JJa <T( , )rW ) 

(ff,T)es„xS„ 

Hence || sdet(S)|| 2 can be written as the sum of the terms 
(4.3.1) 



1 2 , n x 1/2 x n .1/2 

[—{J (sgnai)(sgnri)(sgna 2 )(sgnr2)f JJ^iW^W ) (II ^(0^(0 J x 

\=1 ' \=1 ' 

(^ai(l)n(l) • ••^ 1 (n)T 1 (n) I -^a 2 (l)r 2 (l) ' ' - (n)r 2 (n) ) 

for all 4-tuples (ui, cr 2 , n, r 2 ) G SVi x S n x ,S n x S n . 
Suppose that the sets of indices 

|((7i(z),ti(z)) : i = 1, . . . ,n| and | (cr 2 (i), r 2 (i)) : i = 1, . . . ,n| 

coincide. This is the case if and only if there is a permutation tv E S n such that 
cr 2 (i) = o"i(7r(i)) and r 2 (z) = t\{tt{i)) for z = 1, . . . , n, in which case the term (4.3.1) 
can be written as 

\ 2 n 

i=l 

(^ai(l)n(l) • ' •X ai ( n ) Tl (^ n ), X - 1 ( 7r (i)) Tl ( 7r (i)) • • -Xo-^t,-^))^ (n{n)))- 

Let us fix ai,ri G >S n and let 7r range over S n . Then the expected value of the 
corresponding sum of monomials is 

c(n,n) A 

n! H°«Ti(0n(i)- 

i=i 

If the sets of indices 

:i = l, ... ,n| and | (cr 2 (z), r 2 (z)) : z = 1, . . . , nj 

in (4.3.1) do not coincide, then there is a term X^j, which is in the right hand side 
of the scalar product, but not in the left hand side. The conditional expectation 
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with respect to Xij is because of (4.1) and the linearity of expectation. Hence 
the expectation of the corresponding term is 0. Summarizing, the expected value 
of || sdet B\\ 2 is 

Ec(n, a) t~t , s . 

, ll a «ri(0ri(i) = c(n, fi) per A. 
(<7i,Ti)es„xS„ ' i=i 

□ 

Algorithm 3.5 and Theorem 3.6 ensure the property (1.2.1) of the estimator. 
Theorem 4.3 implies (1.2.3) and (1.2.2) (up to a constant c(n, /z); for sufficiently 
generic measure and scalar product the constant is strictly positive). One can see 
that (1.2.4) is satisfied as well. 

5. A Series of Estimators Conjectured to be Asymptotically Exact 

Let us fix a positive integer d and let A = Mat(d, R) be the algebra of real dx d 
matrices. We introduce the scalar product (•, •) on A by letting 

(a, b) = Tr(a&*). 

Let \x be the standard Gaussian measure on A with the density 

V(a) = (27r)- d2 /2 e -IHI 2 /2. 

Thus to sample a random matrix a G A, we sample each of its d 2 entries indepen- 
dently from the standard normal distribution. 

The construction of Section 4 gives us a permanent estimator. Although we 
don't know the value of c(n, n), we don't really need it, since we can approximate 
it by applying the same algorithm to the identity matrix. 

Hence, for each positive integer d, we obtain a permanent estimator. 

(5.1) The d-th estimator. 

Input: A non-negative n x n real matrix A = (a^). 
Output: A non-negative number a approximating per A. 

Algorithm: Sample n 2 matrices My of size d x d from the standard Gaussian 
distribution in Mat(<i, R). Define nxn matrices B = (6^-), E = (e^) with bij , G 
Mat((i, R) for i, j = 1, . . . , n as follows: 

J i*« if i = i 
Oy — Uijy/aij and — < . . 

10 if i ^ j. 

Apply Algorithm 3.5 to compute d x d matrices sdet 5 and sdet E. Compute 

a = || sdetS|| 2 /|| sdet£|| 2 . 
ll 



Output a. 

We conjecture that the output a approximates per A within an exponential factor 
7^ and that 7^ approaches 1 as d grows. More precisely, we conjecture that there 
exist 

a sequence {7d} of non-negative real numbers such that lim^ — >+oo Id = 1 
and 

a sequence of functions {fd(n, e)} such that for any e > and any d we have 
lim n _^ +00 / d (n, e) = 0; 

so that for any n x n non-negative matrix A and any e > 0, the output a of the 
d-th estimator (5.1) satisfies: 

P {hd-e) n pevA<a< (<y d + e) n per a} > l-f d (n,e). 

Note, that the complexity of the <i-th algorithm is 0(n d2+s ). We discuss some 
plausible reasons why the conjecture might be true. 

(5.2) Why there might be sharp concentration. 

From Theorem 4.3 and Chebyshev's inequality, one deduces that the value of 
c _1 (n, sdetS|| 2 is unlikely to overestimate per A: 

P {c _1 (n,/i)|| sdetSf > K per A} < 1/K 

for every K > 0. Hence the main problem is to prove that c _1 (n, sdet5|| 2 is 
unlikely to underestimate per A as well. 

Suppose that in Algorithm 5.1, instead of the symmetrized determinant of B, 
we take the Cayley determinant Cdet (see Section 1.3). One can prove that 

E || Cdet B || 2 = cT per A 

Moreover, the method of [Barvinok 99] carries over and one can prove that for any 
1 > e > 

P |cT n || Cdet B|| 2 < (ec d ) n peryl) < — 
^ > nln e 

where constant Cd is defined as follows (see [Barvinok 99]): let £1, ... , be inde- 
pendent random variables having the standard Gaussian density (27r) _1 / 2 e -x I 2 . 
Then 



c d = exp{ E ln( ^ + - rf + ^ )}. 



We have lim^ >+00 Cd = 1 and, in fact, Cd = 1 + 0(l/d). 

The symmetrized determinant sdet B can be considered as the average of the 
n! Cayley determinants of the matrices obtained from B by permuting rows in all 
possible ways. It seems quite plausible to the author that the concentration for 
|| sdet B || 2 should be at least as sharp as for || Cdet£?|| 2 . 
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